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Abstract This paper illustrates the application of Lie transform normal-form theory to the construction of 
the 1:2 resonant normal form corresponding to a wide class of natural Hamiltonian systems. We show how 
to compute the bifurcations of the main periodic orbits in a potential with double reflection symmetries. 
The stability analysis of the normal modes and of the periodic orbits in general position allows us to get 
overall informations on the phase-space structure of systems in which this resonance is dominating. As an 
example we apply these results to a class of models useful as galactic potentials. 

PACS. 

45.20.Jj, 47.10.Df Hamiltonian mechanics 
98.62.-g Properties of Galaxies 

1 Introduction 

The 1:2 resonance or 'Fermi resonance' plays a prominent role in nonlinear Hamiltonian dynamics. In galactic dynamics 
it appears in several fashions [1]: to mention a few, in axisymmetric prolate systems it determines the bifurcation of 
the inner thin tube orbits [11]; in triaxial systems with ellipsoidal non-singular cquipotentials it gives the bifurcation 
of banana and anti-banana orbits in the symmetry planes [6,17]. Its interest is clearly not limited to this field and its 
investigation in theoretical and applied nonlinear dynamics has been very active [26,7]: an example is the so-called 
spring-pendulum [9] and an application in satellite attitude dynamics is the tethered system [21]; in chemistry it is 
quite relevant in molecular vibrations [22,23] and in quantum physics for the semi-classical approximation of atomic 
nuclei [4]. 

In [24] we have investigated the relevance of the 1:1 resonance in galactic dynamics in the cases of one and two 
reflection symmetries. Here we want to show how Hamiltonian normal forms can be used to get qualitative and 
quantitative information on the bifurcations connected with the 1:2 resonance. We limit the analysis to systems with 
reflection symmetry with respect to both degrees of freedom: in this case we should more correctly speak of 2:4 
resonance, in view of the structure of the resonant Hamiltonian [12]; however we keep the more standard notion of 
symmetric 1:2 resonance. On this ground, although this case retains all the characteristics of a low-order resonance, it 
also presents some aspects of higher-order resonance systems [2,3] and again makes the results of particular relevance 
for galactic motions. In low-order resonances the bifurcation of new families of periodic orbits is related with loss 
of stability of one of the normal modes, whereas in higher-order resonances the new families appear in a resonance 
manifold from the breakdown of a resonant torus. The peculiarities of the bifurcation sequences from the normal 
modes are worthy of note by themselves, since the approach followed to study systems with a single symmetry like 
the spring-pendulum, is not able in the case of double symmetry to unveil the generic behavior of the system [26]. 

In this work we go a little bit in this direction by considering a generic perturbation up to the degree necessary to 
include resonant terms. Even if this is still not enough to deduce the general behavior for arbitrary perturbations, it 
allows us to gather complete informations on the bifurcation structure near the resonance of the truncated system. In 
fact, we will see that a truncation of the normal form at the first term incorporating the resonance is able to capture the 
essential features of the bifurcation showing how the inclusion of higher order terms in the perturbation is necessary 
to remove some degeneracies. We apply the procedure of the Lie transform normalization whose algorithmic structure 
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helps in the application of higher-order perturbation approaches. Referring to applications in galactic dynamics, where 
several numerical and analytical investigations are available [17,18,19], we briefly present results concerning a class of 
systems with elliptical equipotentials including the well known cored logarithmic potential [6,25]. 

The plane of the paper is as follows: in section 2 we introduce the procedure to construct the approximating 
intcgrable system by recalling the method of the Lie transform [13]; in sections 3 we apply this approach to investigate 
general aspects of the dynamics obtaining second-order estimates of the bifurcation thresholds of the 1:2 periodic 
orbits; in section 4 we analyze the stability of normal modes and periodic orbits in generic position; in section 5 we 
apply these general result to the case of systems with elliptical equipotentials and finally in section 6 we discuss open 
problems and possible further developments. 



2 The model and its normal form 



Suppose the system under investigation is given by a natural Hamiltonian 

H(x,y,p x ,p y ) = -(pl+pl)+V(x,y) (1) 

where V is a smooth potential with an absolute minimum in the origin and symmetric under reflection with respect 
to both coordinate axes. We assume the potential to be expanded as a truncated power series 

N 

V(x,y) = J2v n (x,y) (2) 

n=0 

where V n is a homogeneous polynomial of degree n + 2. The truncation order N is determined by the problem under 
study. In force of the reflection symmetries, the 'zero' order term can be written as 

V a = l -{ulx 2 +Loly 2 ) (3) 

and the odd order terms are all zero. The two coefficients of the quadratic term are written so to represent the 
linearised harmonic frequencies. A system with a potential of type (2) can be treated in a perturbative way as a 
non- linear oscillator system and the two frequencies of the unperturbed system are precisely given by oji and w 2 . In 
order to put the system in a form suitable for a perturbative approach, we perform the scaling [26] 

(x,y,p x ,p y ) -> £~ 1 {x,y,p x ,p y ), e>0 (4) 
and also rescale the Hamiltonian (1) according to 

H = e 2 H. (5) 

Thus, we obtain 



H(x,y 7 p x , Py ) = ^(pl+p 2 y ) + ^(ij 2 1 x 2 +Ljiy 2 ) + J2e 2l V 2j {x 7 y). (6) 

In this way the terms of the expansions are ordered in powers of the small perturbation parameter. 

In general, neither the original system (1) nor its expansion (6) are integrable. However, in several cases the 
dynamics around the equilibrium are regular, namely the measure of chaotic orbits is exponentially small in the 
perturbation and the features are similar to those of an integrable system in a large fraction of phase space. Therefore 
we proceed to construct a 'normal form' [10] for the system, namely a new Hamiltonian series which, in the case of 
2 degrees of freedom, is an integrable approximation of the original one. Its structure is suitable to capture the most 
relevant orbital features of the system. 

In particular, the normal form is 'non-resonant' when the two harmonic frequencies oj\ and ui2 are two real numbers 
so that their ratio is not rational: in this case the normalization produces a 'Birkhoff ' Hamiltonian depending only on 
actions. Therefore, since the new Hamiltonian do not depend on angle variables, we get explicit formulas for actions 
and frequencies of the box orbits parented by the x- and y-axis periodic orbits (the 'normal modes') 

We can instead assemble a 'resonant' normal form by assuming from the start a rational value for the ratio of the 
harmonic frequencies: this assumption produces the presence in the new Hamiltonian of resonant terms, namely terms 
depending on a linear combination of angles with integer coefficients. This expedient is legitimate because, even if the 
unperturbed system is non-resonant for the real value of the frequency ratio 
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P = uj 1 /oj 2 , (7) 

the non- linear interaction between the degrees of freedom, induced by the coupling terms of the perturbation, produces 
a resonant value of p. Its commensurability ratio, say m/n with m, n G N, is determined by the local ratio of oscillations 
in the two degrees of freedom. This in turn is responsible for the birth of new orbit families bifurcating from the normal 
modes or from short-periodic orbits generated by the lower-order resonances. The trick is then to assume that our 
system is such that the ratio (7) is not far from a rational value and then to approximate it by introducing a small 
'detuning' S so that 

p = m/n + S. (8) 

Afterwards we proceed like if the unperturbed harmonic part would be in exact m:n resonance by treating the remaining 
part as a higher order perturbation. We speak of a detuned m:n resonance. 

Let us proceed with the case of the m = 1, n = 2 resonance, so that (8) translates into 

keeping in mind we are in presence of reflection symmetries about both axes: we will shortly see that, in this case, the 
normalization procedure must be pushed at least to the fourth degree in e [5] . To give the system a structure suitable 
to apply the normalization procedure, we perform the transformation 

p x p v 

xi = pi = ——, x 2 = V^y, p 2 = ~t±= (10) 

y/ui V^2 

and introduce the complex variables 



z 1 =p 1 +ix 1 , w 1 =p 1 -ix 1 , 

Z 2 =p 2 +lX 2 , W 2 =p 2 -ix 2 . 

Since we include the detuning term in the perturbation, we introduce a rescaled detuning parameter 5 such that 
6 = Ss 2 . Thus, by redefining the Hamiltonian according to the scaling 

H = 2pH = (1 + 25) ii (12) 
and collecting terms in e, we put the rescaled Hamiltonian into the form 

N/2 

H(z,w:S)=Y,£ 2j H2j(z,w;S) (13) 

3=0 



where the unperturbed term is given by 



and the detuning term 



H (z, w) = ^WiZi + w 2 z 2 (14) 



s 2 S(xf +p() = e 2 5w 1 z 1 , (15) 

rather than being treated as a term of order zero, is considered as a term of order 2. 

The system is now ready for a standard resonant normalization (see [10,27]). Here we are going to apply the method 
based on the Lie transform [10]. An account of the procedure has been given in [24]; we briefly recall here the main 
ideas to adapt it to the symmetric case. The starting point is to treat the canonical transformation as a 'flow' along 
the Hamiltonian vector field associated to a generating function. Let us consider a phase-space function expanded as 
a power series in the canonical variables 

N/2 

G = J2? 2l G2i (16) 

i=i 

where the odd-order terms are zero in force of the reflection symmetries. Since we want the method to work with the 
series expansion given in (13), we assume that the non- vanishing terms are polynomials of the form 
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G 2 i(z,w;S) = ^2d J g 2 (i- j+ i)(z,w), g m eV m , to = 0, ... Z. (17) 
j=o 

where P m denotes the space of homogenous polynomials of degree to + 2 in the (z, w) coordinates which are invariant 
with respect to the reflections 

{z 1 ,Z 2 ,W 1 ,W 2 ) -> (--Zi, Z2,-«>1,W2) (18) 

(z 1 ,z 2 ,w 1 ,w 2 ) -> (zi,-z 2 ,wi,-W2) (19) 
and their combinations. To the series (16) is naturally associated the linear differential operator 

e£G =E^G, (20) 

k 

whose action on a generic function F is given by the Poisson bracket: 

£ G F = {F,G}. (21) 

The original Hamiltonian system (13) undergoes a canonical transformation to new variables (Z,W), such that 
the new Hamiltonian is 

K(Z,W) =c Cg H(z,w), (22) 
where K is assumed in the form of a series expansion similar to (13), namely 

AT/2 

K(Z,W)=Y,^K 2J (Z,W). (23) 
j=o 

To construct K starting from (13) is a recursive procedure exploiting an algorithm based on the Lie transform [10, 13]. 
To understand how it works, let us consider the first step of the procedure, namely let us perform a transformation 
given by a function G 2 considered as the first term in the generating function (16). The general relation (22) takes 
the form 

K + e 2 K 2 + ... = (1 + C G2 + ...)(H + e 2 H 2 + ...). (24) 
By equating polynomials of the same order in e, we get the system: 

K = H (25) 
K 2 = H 2 -C G2 H Q (26) 

K 4 = H 4 -C G2 H 2 -lc% 2 H (27) 



K n — C G2 H + R n (28) 

where the 'rest' R n contains terms which are known if the preceding n — 1 equations have been solved. 

Equation (25) simply states that the zero order new Hamiltonian coincides with the zero order (unperturbed) one. 
To proceed, we have to solve the second equation to find K 2 , a differential equation involving two unknown functions, 
K 2 and G 2 . To overcome this problem we have to make some decision about the structure the new Hamiltonian K 
must have, that is we have to choose a normal form for it. We make the choice that K has to satisfy 

{K,H } = (29) 

or equivalently 

C Ho K = 0. (30) 

In this way the system with Hamiltonian K admits H as a new integral of motion. To satisfy condition (30), let us 
consider for greater generality, the n-th step in the normalization procedure and impose that K n and G„ be solutions 
of the system 
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We rewrite the first equation as 



C Ho K n =0. W 



£H G n + K n = R n . (32) 



This is the so called homological equation. By the way, in force of the reflection symmetry, for odd order terms equation 
(32) gives the trivial solution K n = G n = 0, proving the consistency of the assumptions concerning the expansions 
(16) and (23). For n = 21, if we recall that the detuning parameter is assumed of order two, equation (32) becomes 

l i i 

3=0 J=0 j=0 

with g m , k m and r m G V m , m = 0, . . . I. Since Ch ■ V m — > V m , to solve equation (33) is equivalent to solving I = n/2 
equations of the type 

72 

^H 92(i-j+i) + k 2 (i-j+i) = r 2 (i-j+i), 3=0,...,-. (34) 

Now, thanks to the semisimple character of the linear operator Ch , its kernel and its range are in direct sum over 
the space V m . This implies that the system of equations (34) and hence the n-th homological equation, can always be 
solved if K n satisfies (30) . 

We have so far constructed the normal form Kq + • • • + K n and the generating function G 2 + • • • + G n : we use 
them to compute K n+2 and G„+2 and so for. However, as we have observed above, the series are divergent, thus we 
must truncate the procedure at some finite order, say M. In the case of a m = 1, n = 2 resonance in the presence 
of reflection symmetries about both axes, the normalization procedure must be pushed at least to order M = 4. 
Therefore, generalizing the assumptions made in [3] , we assume that the non- vanishing terms in the series expansion 
of the original Hamiltonian are given by 

H = l -{x\+p\) + {xl+pl) (35) 

H 2 = 5(x\ + p\) - ^{ax\ + bx\ + cxlxl) (36) 

Ha = - (aixf + b\x\x\ + c\x\x\ + dix®) , (37) 
6 

where the arbitrary coefficients appearing in the higher order terms represent the most general potential truncated to 
degree six in the coordinates and complying with the enforced double reflection symmetry. The choice of coefficients 
and signs is suggested by the values the coefficients take in the commoner physical cases: we remark that, in applying 
the results obtained in the following to specific model problems, we have to take into account that the Hamiltonian 
(35-37) is in the form 'prepared' for normalization. Therefore the canonical variables are rescaled according to (10) 
and the frequency ratio is expanded in series of the detuning as in (9). These transformations affect the numerical 
values of the coefficients of the various terms. 

The outcome of the normalization procedure described above is a series of the form (23) that, by introducing 
action-angle like variables by means of 



Z = iV2J c 
W = -iV2J e 



(38) 



can be written as 



K(J,8)=J2e 23 K 23 (J,8), (39) 

3=0 

which, coherently with the truncation order M = 4, implies that a reminder of order six in e is neglected. In these 
canonical variables it is possible to express the series in a manageable form whereas the use of Cartesian variables 
(both in real or complex form) gives quite cumbersome formulas. The non-vanishing terms in the normal form turn 
out to be the following: 
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K = J 1 + 2 J 2 (40) 
1 

-< 
4 



K 2 = 2~6J 1 -laJ 2 1 -]cJ 1 J 2 -hjl (41) 



= {-f/ + ^ ai ) Jl + {-^-/ + J3 + _L (_ 186c _ ^ + 48ci ) J x J 2 2 

+ jj 2 J 2 [- (36a + 5c) c + 48&i + (3(a - c)c + 86i) 008(461! - 20 2 )] . (42) 
192 

The resulting Hamiltonian (39) is the basis for all subsequent work. We observe that up to order 2, the dependence is 
only on action-like variables: therefore, up to this order they are true conserved actions with angles evolving linearly 
on the invariant tori. The resonant combination of the angles with zero phase difference appears only in the term of 
order 4 (and higher) and this is the reason why this is the minimum order required in this study. In ref. [3] it is proven 
that this phase value always occur in potential problems and that in more general Hamiltonian systems the phase can 
be different from 0, ±7r. 

The resonant terms determine a nonlinear coupling of the two degrees of freedom. Although the system is integrable 
by construction, the solution of the canonical equations is in general not expressible in terms of elementary functions. 
We can however exploit the second conserved quantity to reduce the system and we show in the subsequent sections 
how to use this technique to gain understanding of the structure of the phase-space. 



3 Bifurcation analysis of the 1:2 symmetric resonance 

The essential information we need concerns the existence and stability of the periodic orbits associated to the resonance 
[5,24]. We can compute the thresholds for the bifurcations sequences in terms of the parameters relying on the regular 
nature of the dynamics given by the normal form. In two degrees of freedom, if a Hamiltonian possesses a second 
independent integral of motion, the system is Liouville integrable. Due to the normalization procedure, we have 
obtained the Hamiltonian (39) with the second independent integral of motion 



K = J 1 + 2 J 2 



(43) 



We can use this integral to reduce the dimension of the problem by performing the canonical transformation to 'adapted 
resonance coordinates' [27] 



( Ji = £ + 2TZ 

J 2 = 2£-K 
ip= 40i - 20 2 
IX= 201+402 



(44) 



It can be easily seen that \ is cyclic and its conjugate momentum is proportional to the additional integral of motion, 
namely 



£ = ( Ji + 2 J 2 )/5. 



(45) 



The value of £ plays a special role in the discussion of the orbit structure and we will refer to it as the distinguished 
parameter [9] . Thus we introduce the 'reduced Hamiltonian' 



K{Tl, i/>;£) = 5£ + e 2 K, 1 {K) +e 4 /C 2 (ft, i/>). 



(46) 



In the computation of (46) and through the rest of the paper the use of algebraic manipulators is almost indispensable. 
By using mathematica ® we obtain the following expressions for JC\ and /C 2 : 



/c 2 (^,v) = 



3b(-2£ + Tl) 2 + 3a(£ + 2TZ) 2 + 2(£ + 2K){2c£ - cTZ 



17 

64 



(a 2 +4b 2 )-^c(a + b) 



1 



10 



(47) 



32 



12 



«i + + ci + — di £ A + A(R) + B(K) cos ip (48) 



with 
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T?3 

A(K) = — (-816a 2 + 516 2 + 288ac - 726c + 20c 2 + 1280ai - 3846i + 192ci 
384 

- 160di) + — (-204a 2 - 51b 2 - 48ac + 426c + 5c 2 + 320 ai + 646 x 

£ 2 1Z 

- 112ci + 160di) - — (306a 2 - 3066 2 + 252ac + 726c + 55c 2 

- 480ai -336&i - 192ci + 960di) , (49) 

B(TZ) = (2£ - n) (£ + 2TZ) 2 (3ac - 3c 2 + 8&i) . (50) 
192 

Considering the dynamics at a fixed values of £, we have that JC defines a one degree of freedom system with the 
following equations of motion 

ip=^ (6b£ - 3c£ - 3bTZ + 4cTZ - 6a(£ + 211) + 165) e 2 + + |j| cos ipj e 4 (51) 

K= ^-((2£- ll)(£ + 2K) 2 sin iP)e 4 (52) 

where 

/i = 3ac- 3c 2 + 8&i. (53) 

The fixed points of this system give the periodic orbits of the original system. 

The pair of fixed points with 1Z = 2£, 1Z = —£/2 correspond to the normal modes, that is to the periodic orbit 
along the x-axis (J 2 = 0) and to the periodic orbit along the y-axis (Ji = 0) respectively. Additional periodic orbits 
may appear when the system passes through the resonance. These periodic orbits 'in general position' exist only above 
a given threshold in the distinguished parameter £ when the axial orbits change their stability. This phenomenon 
can be seen as a bifurcation of the new family from the normal mode when it enters in 1:2 resonance with a normal 
perturbation. The phase between the two oscillations determines the nature of the families: they are respectively given 
by the conditions ip = (banana orbits) and ip = ±ir (anti-banana orbits), where we use the nicknames introduced in 
ref.[17]. These phase conditions are solutions of 1Z = (when 1Z ^ 2£ and 1Z ^ —£/2) and determine the corresponding 
solutions of ip = 0. 

Let us start looking for banana orbits. Following a standard approach [14], we set tjj = in the equation ip = 
and look for a solution in the form 

K = K +K l s 2 + 0(8^) (54) 

so that the righthand side of (51) vanish up to fourth order in e. For ip — 0, we substitute (54) in (51) and collect 
terms of the same order in e. Equating to zero the coefficient of the second order term in e, we find that TZq has to 
satisfy 

3(26 - 2a - c)£ + (4c - 36 - 12a)"K + 165 = (55) 
This equation admits solution only if 



v = 12a + 36 - 4c ^ (56) 



and, if this condition is satisfied, we find 



^ -3(2a-2b + c)£ + 165 
n *- n z°- 12a + 36 -4c ' (57) 

Once IZq is computed, we look for IZi such that the fourth order term of the righthand side of (51) vanish. Since (56) 
is satisfied, we find a solution 

H 1 =Hbi(£;6). (58) 

The corresponding fixed point is given by 



K = K B =n B o+KBie 2 , V = 



(59) 
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and determines the banana orbits (there are two of them): 



J 1 = J 1B =£ + 2TZ B (60) 
J 2 = J 2B = 2£ -K B - (61) 



Similarly, if (56) is satisfied, for 4#i — 29 2 = ±7r, we find 



,h = J 1A =£ + 2K A (62) 
J 2 = J 2A = 2£ -TZ A , (63) 

which correspond to the anti-banana orbits. In view of (44), the constraints 

5£ 

0<Ji<5£, 0<J 2 <y, (64) 

applied to these solutions give the condition of existence for these periodic orbits in general position. Whether these 
conditions are satisfied or not, it depends on the parameters of the system. For the banana orbits we find that at the 
zero perturbative order 

_ 5(36 - 2c)£ + 325 
13 ~ 12a + 36 -4c ' 
5(6a — c)£ — 16S 

J ™ = 12a + 36-4c • (65) 

Thus, we get different existence conditions according to the sign of the constant v defined in (56). Namely, taking e 
small enough such that the constraints (64) remains satisfied up to the first perturbative order, banana orbits bifurcate 
in the following cases: 

* & <°^*<»-*> ::£Sife <"<*„■ < ct > 

where the critical values 



16 - 16 (306a 2 - 33ac - 8c 2 - 480ai + 56&i) ~ 9 , , ^ 
£ ™ = Wa^) 5 - 15(6^ ^ + 0(e*) (68) 

. 32 - 16 (153b 2 -72bc- 20c 2 + 192^-480^) - 4 
f 52 = 5(^36) ^ 15(36 -2c)3 * e + 0(e ) (69) 

correspond to the solutions of J 2 b — and Jib = and respectively determine the bifurcation of the banana orbits 
from the .t— normal mode and the y— normal mode in the first sub-cases of (66) and (67). In the other two sub-cases 
the periodic orbit bifurcating from one of the normal modes disappears on the other. 

A similar argument provides the existence condition of anti-banana orbits. Since to the first perturbative order 
J\A = Jib & n d J 2 a = J 2 b, the birth of anti-bananas is given by the same conditions on the coefficients given above. 
However the higher order terms in JkA, JkB, fc = 1, 2 are in general different: the discrimination between the thresholds 
for the bifurcation of the two families is possible only by going at least to second order [5]. This result was expected 
on the basis of the structure of the new Hamiltonian (39). 

Nevertheless, up to the second order in e we find that the threshold value corresponding to the bifurcation of anti- 
banana orbits from the y— axis, £a 2 , coincides with the critical value £b 2 of (69). Hence, if banana and anti-banana 
orbits bifurcate from the y— normal mode, they do it concurrently. On the other hand the bifurcation from the a;— axis 
orbit occurs at 

£=£au 5(6a-c)>0, (70) 



where 
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c . 16 - 16 (306a 2 - 39ac - 2c 2 - 480 ai + 406i) ~ 2 2 : 

£A1 = 5(6^) ^ 15(6o-c)3 ^ £ + °( £ ) 

which, at second order, is different from (68). Thus, we obtain the following existence conditions for anti-banana orbits: 



^-0and^- 2c) <0,{;;| ::: )^;^^ <fAi , (73) 

with H* 1 ^* 1 ' *jj">° (74) 
{£ A1 < £ m , if 5fi< 

We can write explicitly the relative magnitude of the two thresholds for bifurcation from the normal modes because 
we find 

£ A 2-£B2 = + O(e 4 ). (76) 

The first of these expressions points out that the hierarchy of bifurcations from the x— normal mode is determined by 
the sign of the constant fi defined in (53). The second confirms that, at the order of our perturbative treatment, the 
bifurcations from the y— normal mode occur simultaneously. Which of the two scenarios is actually happening depend 
on the value of the parameters according to the conditions listed above. 

In Fig.l we plot the surfaces of sections relevant for the first scenario concerning the bifurcations from the x— normal 
mode: the upper panels are the (x,p x ) sections after the bifurcation of the stable banana (left) and after the subsequent 
bifurcation of the unstable anti-banana (right); the lower panels are the corresponding (y,p y ) sections. This occurrence 
is typical of a lower order bifurcation: the normal mode losing stability at the first bifurcation is the origin in the 
upper panels and is the last contour in the lower. The contour in the lower left panel is the unstable mode that regains 
stability in the right panel. 

The second scenario in which the two families appear together from the y— normal mode is the typical occurrence 
in higher-order resonances: pictorially it correspond to a transition from the non-resonant invariant tori around the 
normal mode directly to the appearance of the resonance manifold as it emerges in the two lower panels with stable 
and unstable orbits arising in pair. 

For practical purposes, e.g. to compare these findings with the outcomes of numerical simulations, it would be 
more useful to have the expression of the bifurcation curves in terms of 'physical' parameters. The most natural way 
to represent the thresholds is that plotting curves in the (E, 5)— plane, where E is the physical energy of the system 
defined by (1) 

H(x, y,Px,P y ) = E (77) 

and 5 is the 'true' detuning defined in (8). According to the rescaling (12), on the x— axis orbit (J 2 = 0, J\ = 5£), we 
have 

e 2 K = 5£e 2 + (- + 10£^j e 4 + • • • = (1 + 2~8e 2 )E. (78) 

The dots are present to recall that a reminder has been neglected. The series from equation (78) can be used to express 
the physical energy E in terms of £ [5,16], namely 

75 

E = h£e 2 - ^a£ 2 e 4 + 0(s 6 ). (79) 
8 

Thus, up to the second order in e, for £ satisfying equations (68) and (71) and 5 as in (9), we obtain the following 
threshold values 

= _U^ S _ 16(41 4a 2 -51 fl; :-8c 2 -48Q fll + 56^ 2 + 

6a — c 3(6a — c) A 

16 16 (414a 2 - 57ac - 2c 2 - 480ai +406i) , , 

Eai = ~, S ^ — - 3 lL{? + (6 3 ), 81 

6a - c 3(6a — c) A 
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Figure 1. Dynamics of the 2:4 doubly-symmetric resonant Hamiltonian, see text for explanation. 



for the bifurcation of respectively banana and anti-banana orbits from the a;— axis orbits. 

A similar argument gives the threshold value for the bifurcations from the y— axis orbit. Since £a2 = &B2, we have 
Ea2 = Eb2 and by using the relation between the true energy and the distinguished parameter on the normal mode 
we get 



E 



B2 



32 16 (996 2 + 366c - 68c 2 + 192ci - 480di) 



2c- 36 



3(36 - 2c) 3 



5 2 + 0(5 3 ). 



(82) 
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4 Stability analysis of the 1:2 symmetric resonance 

Let us now consider the question of the stability of periodic orbits: this analysis complements that of the previous section 
allowing us to test the relation between change of nature of normal modes and bifurcation of a new family. For banana 
and anti-banana orbits an ordinary investigation of the equations of variations of the system is enough to perform the 
linear stability analysis. However, in the case of axial orbits, action angle variables have singularities on them and this 
also affects the adapted resonance coordinates. However the remedy to this problem is quite straightforward: to use a 
mixed combination of action angle variables on the orbit itself and Cartesian variables for the other dof. 

Let us start with the stability analysis of the periodic orbits in general position. We have to investigate the fate 
of a normal perturbation of the periodic orbit under test. The system of differential equations for the perturbations 
(dtp, 51Z) is given by 

= K r n ){lt)- ( 83 ) 

dt J \—K>w -K-iiip J \oTZ J 

Here we again use the reduced Hamiltonian (46) and, with a small abuse of notation, we assume without denoting it 
explicitly that the entries in the Hessian matrix are evaluated on the specific orbit we are interested in. Then, the sign 
of the determinant 



A(K,fr£,S) = K^+-KwilCw (84) 

computed on the periodic orbit determines the fate of the perturbation: if A(1Z, tp; £, 5) is negative it gives the frequency 
of bounded oscillating solutions thus determining stability; a change of sign, as a consequence of varying £, produces 
a stability transition. 

On the banana and anti-banana orbits we respectively have 



Alm „. A . (30a£ -5c£ -I65)(l5b£ -10c£ + 32S) 2 * Rn 

A(H B ,0) = A B = ^ 768(12a %_ 4c)2 £ + <>&) (85) 

(30a£ - 5c£ - 16d)(l5b£ - 10c£ + 32S) 2 fi nl 8 , 
A(K A ,,) = A A = 768(12a + 36-4 c)2 ~ £ + ° {£) (86) 

and thus we see that the parameter fi plays an important role also for stability. Comparing with (68) and (69), for 
dfi > banana orbits are stable in the case they bifurcate from the x-axis orbit (£ > £bi, if Sv > and SN > 0) 
and unstable in case their bifurcation occurs from the y-axis orbit. Otherwise, we have instability (stability) when the 
bifurcation occurs from the x- normal mode (j/-normal mode). Since A a = —As up to the third pcrturbative order, 
anti-banana orbits turn out to be unstable when banana orbits are stable and viceversa. Actually, the fourth order 
terms in (85) and (86) are different, but their difference is again a multiple of fj,. 

We have also seen in (75) that the bifurcation order from the x-axis depends on the sign of Sfi. Thus, we can now 
state that, if Sfi is positive, we have at first the bifurcation of (stable) banana orbits followed by (unstable) anti-banana 
orbits. On the contrary, for negative values of 5/j, the bifurcation order and stability nature are inverted. 

Let us now study the stability of the normal modes. Considering the a;— axis orbit, we use action-angle variables 
on the orbit and Cartesian variables on the normal bundle to it, namely 



(87) 

[ Py = V 
so that the periodic orbit is given by 

Y = V = 0, J = 5£. (88) 
In these coordinates, the system of differential equation for the perturbations of the normal mode is given by 




- ( 5Y \ - ( K vy Klvv \ ( SY 
dt \sv) - \-K YY -iCyvJ \sv 



(89) 



where K, = K.{Y, V, 9, J). However the matrix of the second derivatives of K, on the periodic orbit depends on 0(f) = ut, 
where 
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dIC (Vba£ ~\ 2 /1275 , 15 



To remove the dependence on time we introduce complex coordinates 



and perform the 'rotation' 



z= Y + iV 
w = Y-iV 



z= Ze- 2iut , 
w = We 2iut . 



In this way, the equations of variation (89) on the periodic orbit (88) give 



±{SZ\_ (An A 12 \ ( 8Z\ 
dt \SW) ~ \A 2 i A 22 ) \5W ) 



where, up to the second perturbative order, 



(91) 



(92) 



(93) 



1 - 25 

An = -A 22 = -(30a£ - hc£ - 16S)e 2 + — (306a 2 - 36ac - 5c 2 - 480 ai + 486i) £ 2 e 4 (94) 

25 

A12 = -A 2 i = - jgg/^V. (95) 

By solving det A = we find, as expected, that the critical values of £ which determines a change in the stability of 
the a:— axis orbit are precisely given by the bifurcation values 

£ = £bi and £ = £ai (96) 

as defined by (68) and (71). Regardless of the the nature of the occurring bifurcations (this is given by the sign of /it), 
the first one produces a transition from stability to instability of the x— normal mode and the second one a return to 
stability. 

Let us now consider the stability of the y— axis orbit. Since the periodic orbits in general position bifurcate 
concurrently from this normal mode, we expect that the y— axis orbit remains stable after the bifurcation. To verify 
this assert, we proceed as above by introducing the coordinates 



(97) 

[Px= U 
so that the periodic orbit is given by 

X = U = 0, J =\ £ - ( 98 ) 
The system of differential equation for the perturbation of the normal mode is given by 

d_ (5X\ _ ( t ux Kuu \ (5X\ 

dt{8u)-{-jc xx -ic xu ){su) W 

where now K = K,(X,U,9, J). Since we are dealing with a perturbation of a 1:2 symmetric resonance, the terms 
proportional to cos(4#! — 29 2 ) in K± are of second degree in Ji and, as a consequence, the matrix of the second 
derivative of K computed on the y— normal mode does not depend on 9. Thus, we do not need to perform the 
transformation (91). The equations of variation (99) give 




d_ ( SX\ _ (fin fli 2 \ ( SX 

dt \su ) ~ I n 2X q 22 ) I su 



(100) 



where 
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fin = fi 22 = (101) 



fii2 = -fi 2 i = 1 - ( — £ + 26 ) - — £ 2 (c(18& + 5c) - 48ci) e 4 . (102) 



Thus, 



'4 / 384 



— (-450bc£ 2 + 25c 2 £ 2 
584 v 

- 960(£ (5 + 15365 2 + 1200£ 2 ci) e 4 + (9(e 6 ) (103) 
has a positive zero order term, which implies that for e small enough the y— axis orbit is always stable. 

5 Application: 1:2 resonance in systems with elliptical equipotentials 

We now illustrate how to apply the above theory to some cases relevant for galactic dynamics. In particular, we are 
interested in the question of the stability of axial orbits in triaxial potentials and of the possible existence of additional 
stable families of periodic orbits. The present analysis of 2 DOF systems is a first step in this program because in 
allows us to study the dynamics in principal planes of triaxial systems with reflection symmetries. In turn these studies 
are useful to solve problems like the construction of self-consistent equilibria, the computation of isophotal shapes and 
velocity ellipsoids, etc. 

We consider a fairly general class of potentials with self-similar elliptical equipotcntial and unit 'core' radius of the 
form [6,11,24] 

!/ 2\ Q / 2 
i(l + .T 2 + ^) , 0<a<2 
ilog(l+x 2 + ^J , a = 0. 

The physical parameters are q, the ellipticity of the equipotentials and a, a shape parameter: on the two extremes of 
its range, a = corresponds to the standard logarithmic potential and a = 2 to the anisotropic harmonic oscillator. 
The family of potentials (104) admits a series expansion of the form (2) with 



2 V q 2 

a - 2 / 2 y< 



V2 = ^[x 2 + y - 2 ) , (106) 



Since the unperturbed frequencies now are uj\ = 1 and u>2 — l/<7, we introduce the detuning parameter (9) 

5 = q-\. (108) 

Performing the scalings (10) and (12) and collecting terms of the same order in e we obtain the Hamiltonian function 
(13), where the non- vanishing terms are now given by 

H = l -(xl+pl) + {xl+pl) (109) 

H 2 = 5{x\ + p\)- 2 -^(x\ + 2xlf (110) 

H 4 = - J(2 - a){x\ - 4x1) + (2 = ^ a) (x\ + 2x\f . (Ill) 
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Hence, in this particular case 



a=i±^(2-a), (112) 

6 = 2(2-a)(l-25), (113) 

c = 2(2-a), (114) 

oi = i(2-a)(4-a), (115) 

6i = 6ai, (116) 

ci = 12ai, (117) 

di = 8ai. (118) 



The presence of terms proportional to 5 in H4 is due to the dependence of the potential (104) on q. However the 
existence and stability analysis of the periodic orbits of the system follows exactly the same way of the preceding 
section. 

Since we have 12a + 36 — 4c = 4(2 — a) > 0, the non degeneracy condition (56) is satisfied and the system is able 
to exhibit the bifurcation of periodic orbits in general position. In order to establish their existence, we look at the 
sign of 6a — c and/or 36 — 2c. We find 



6a-c=(l + 6J)(2-a), (119) 

thus for 5 > these systems fall in the cases (66) and (72). Hence banana and anti-banana both bifurcate from the 
x— axis orbit respectively for E = Ebi and E = Eai, where the thresholds in terms of the physical parameters are 
given in (80,81) and for potential (104) turn out to be 



2 



16 / 1\ 8(41a - 10) / 1\ 



2 



16 / 1\ 8(53a + 14) / 1\ 
^ 1= 2^-2j + i2^(*-2> • (121) 



Since in this case we have 



^ = 3(4-cv 2 ) > 0, (122) 
in agreement with the general expression (75), the difference between the two thresholds is 

^i-^i=32^L(g-i) a . (123 ) 

This relation and equations (85) and (86) establish the bifurcation of stable banana orbits followed by unstable anti- 
banana orbits. These results generalize those already obtained in the work on the logarithmic potential [5,16] and 
provide good approximations to the numerical investigations available in literature [17,18]. 



6 Conclusions 

We have presented the investigation of a fairly large class of natural reversible Hamiltonian systems close to the 1:2 
resonance and endowed with reflection symmetry with respect to the configuration variables of both DOFs. By means 
of a resonant detuned normal form, we have obtained a general description of the bifurcation scenario of periodic 
orbits in general position (banana and anti-banana) from the normal modes. 

Equations (80-82) provide the energy levels for these bifurcation in terms of the detuning and the other physical 
parameters characterizing the system. We have found that the coefficients of the quartic term in the potential essentially 
determine the distinction between the bifurcations from cither normal mode: our main result is that if banana and 
anti-banana orbits bifurcate from the y— normal mode, they do it at the same energy level; if instead the bifurcations 
are from the x— normal mode, their sequence (and the stability of the new orbits) is determined by the parameter /j, 
defined in (53). We remark that for the reliability of the predictions based on these parameters it is essential that in 
the computation of the normal form terms of order four (degree six) in the potential are included. As an example, in 
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the case of the 'a— models' (104), if we arrest the expansion of the potential at order two and normalize up to order 
four, we get 

H= -9 (4 -a 2 ) < 0, (124) 

obtaining a wrong prediction for the bifurcation sequence. 

However, although these results are in excellent agreement with numerical simulations, we can not deduce from 
them the generic behavior of the system. For example, it is difficult to say if the concurrent bifurcation from the 
y— normal mode is persistent and, if not, at which order it split. Actually, the 'catastrophe germ' for a 1:2 symmetric 
resonance is given by the second order term in the perturbation, but has infinite codimension [8]. In simpler terms, 
this means that truncating at order four could not be enough and 'a priori' one has to add (infinitely) many higher 
order terms to the series expansion (13) to obtain a faithful description of the true dynamics of the system. Additional 
efforts are therefore necessary for a full understanding of this problem. 
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